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Abstract 

We consider theoretically a current flowing perpendicular to interfaces 
of a spin-valve type ferromagnetic metallic junction. For the first time an 
effective approach is investigated to calculate a simultaneous action of the 
two current effects, namely, the nonequilibrium longitudinal spin injection 
and the transversal spin-transfer surface torque. Dispersion relation for 
fluctuations is derived and solved. Nonlinear problem is solved about 
steady state arising due to instability for a thick enough free layer. 



1 Introduction 

Great attention is attracted now to features of current flowing through ferro- 
magnetic junctions of a spin-valve type, i.e. structures with contacting ferro- 
magnetic thin layers, one of them having pinned and the other free spins. As 
experiments showed, current can influence substantially the free layer magnetic 
state of such junctions that leads to resistance jumps pQ- as well as mi- 
crowave emission 0]- [H]. 

Mechanisms of the current effect are not completely clear so far. A mech- 
anism was proposed in 7 of the current effect on the free ferromagnetic layer 
magnetization M due to injection of nonequilibrium longitudinal (i.e., collinear 
to M) spins into the layer. Detailed theory of the mechanism was developed 
in |10j . The injection creates a nonequilibrium carrier spin polarization in 
the layer. This polarization, in its turn, contributes to the sd exchange en- 
ergy U s d(j) an d to the corresponding sd exchange effective field H sc j(j), which 
is dependent on the electric current density j. For current densities exceed- 
ing some threshold, a reorientation first-order phase transition occurs and the 
magnetization vector M direction may change abruptly. Such a current induced 

*e-mail:zil(Q! ms.ire.rssi.ru 



1 



magnetization reversal (or current induced switching) leads to resistance jumps, 
which corresponds to the experimental data p^- [3]. 

On the other hand, another mechanism of current effect on the ferromag- 
netic layer state was proposed ^fl El l° n g before appearing of the experimen- 
tal results P0~[ni- According to the mechanism, the transverse (with respect to 
M) mobile electrons spin current is to vanish near the interface between two 
(pinned and free) contacting non-collinear ferromagnets. It occurs due to sta- 
tistical velocity spread of the mobile electrons not because of the relaxation. 
The electrons interact with magnetic lattice (sd exchange interaction) and that 
is why a torque appears at the interface. This spin torque acts on the lattice 
and transfers the lost spin current to it. This spin-transfer torque may exceed 
the dissipation torque at large enough current densities. Then the initial state 
becomes unstable and the magnetization reversal occurs. As estimations show, 
such a switching mechanism can agree with experimental results PP-® also. 

In this connection, a very significant question arises: how are things going in 
the real experiments when both effects coexist, namely, longitudinal spin injec- 
tion and corresponding current dependent effective field, and current dependent 
spin-transfer torque at the interface of the magnetic junction. Up to now, these 
effects were studied separately. Meanwhile, these effects do not only coexist, but 
influence each other also. Therefore, both effects are to be taken into account 
simultaneously, in scope of a unified theory, to understand better the experi- 
mental situation. Such a theory was developed for the first time in our previous 
preprint an d article 14 . The aim of the present paper is to elaborate a 
new effective approach for solving the problem. This approach allows not only 
get all the previously achievable results (see p3E]), but additionally to solve 
nonlinear problem about the steady state arising due to instability development. 
The approach based on the interesting feature of the effective field H s d that was 
revealed in |13M14| . This field is created due to the injection of nonequilibrium 
longitudinal mobile spins into the free layer volume and depends on the bulk 
parameters of the free layer. But a point of the field localization occurs near 
the interface. This feature suggests a possibility to include an effect of the field 
into boundary condition and shows the transversal spin flux becomes, in fact, 
discontinues at the interface. The simplification arises because of no singularity 
remains in the equations of motion. 

2 Model 

We will consider a spin- valve type magnetic junction with current flowing across 
the layer interfaces (CPP geometry). Ferromagnetic metal layer 1 has pinned 
orientations of the lattice and mobile electron spins. Another ferromagnetic 
metal layer 2 contacts with the layer 1 at a point x = and has free lattice 1 
magnetization M and mobile electron magnetization m, so that the magnetiza- 
tions direction can be changed by an external magnetic field H or spin-polarized 
current density j. There is a very thin nonmagnetic spacer between the layers 

1 We mean a lattice of magnetic ions considered in a continuum media approximation. 
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1 and 2, which will be considered as a geometrical plane x = . To close the 
electric circuit, a nonmagnetic metal layer 3 exists in the region x > L. 

External magnetic field H, lattice and mobile electron magnetizations of the 
layer 1, Mi and mi, as well as corresponding magnetizations of the layer 2, M 
and m , are assumed lying in the junction plane x = Q. Vectors Mi and M 
make an angle \ of an arbitrary value under current j = 0. After the current is 
turned on electrons transfer from the layer 1 into the layer 2, that is j/e > 0. 
Then the electrons appear in a non-stationary quantum state in the layer 2 and 
"walk" between the spin subbands of this layer. It corresponds to precession of 
the mobile electron magnetization around the lattice magnetization M. We may 
introduce two components of the vector m by an equality m = m|| +mj_, where 
longitudinal component m|| is parallel to M and transversal component mj_ is 
perpendicular to M. As it was first shown in (lit I12j . angle of the precession 
decreases with coordinate x increasing and tends to zero at x w A_f, where 
A_f ~ 1 nm is an electron quantum wave length at the Fermi surface. This may 
be valid only because of electron velocity statistical spread and no relaxation 
processes are needed to provide such a behavior (more details may be found in 
the recent preprint also). Therefore, both components m|| and m_L exist in 
the region < x < Xp, but only component my remains in the region x > Xp 
(m ± =0). 

The region < x < Xp was introduced firstly by Slonczewski JT] and Berger 
[T2| and further we will refer to it as " SB layer" . We follow the original works 

IT2] and assume ballistic regime of motion into the SB layer. The validity 
criteria for this assumption may be written as Xf < l p , where momentum free 
pass length l p may be estimated as l p ~1-10 nm for metals at room temperature. 



3 Equations 

We describe the motion of vector M by means of the Landau-Lifshitz-Gilbcrt 
(LLG) equation [T5] : 



-^r+7[M,H e// ]-- 



M 



dM 



0, 



(1) 



where 7 is the gyromagnctic ratio, k is dimcnsionlcss damping constant (0 < 

K < 1), 

,a 2 M 



H 



eff 



H + H Q 



A- 



dx 2 



H d + H 



(2) 



is effective field, H a = /3M is anisotropy field, A is the intralattice inhomo- 
geneous exchange constant, is demagnetization field, H sc i is sd exchange 
effective field. The latter takes the form ITS1 



H s d{x,t) 



SUsd 
~SM{x,t)' 



(3) 
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where S(. . .)/8M.(x,t) is a variational derivative, and U s d is sd exchange energy, 



U sd = -a[ m(x',t)M(x',t)dx', (4) 
Jo 

parameter a being a dimensionless sd exchange constant (typical estimation is 
a ~ 10 4 -10 6 Due to the last term in (|2J), the motions of M and m vectors 
appear to be coupled. 

We describe the motion of vector m by means of continuity equation (e.g., 
see PSE]): 

5m dJ r Am 

— + - +7 a[m,M] + — = 0, (5) 

where r is a time of relaxation to the local equilibrium value m = mM, Am = 
m m = AmM and M = M/M is the unit vector, J is mobile electron 
magnetization flux density. 



4 Vector boundary conditions 

The SB layer may influence as a special boundary condition if the thickness 
of the layer 2 is large enough, namely, L \p. Moreover, the value of Xp 
should be assumed the smallest one among the other lengths in the system 
(e.g., \f~A ^> \p and I 3> Xf)- We may sum the Eqs. and © and so 
consider the continuity equation for the total magnetization vector M + m. Let 
us integrate the sum equation over the region < x < Xp , the integration being 
denoted as 

/Af+e 
(...)dx = {...). (6) 

Then the following three types of summands appear. The first type sum- 
mands are proportional to small length Xf and may be neglected for actual 
values of parameters (fields: H, H a , H^; relaxation parameters: k, t; and fre- 
quencies w in a microwave region or lower) . The second type of summands may 
be presented as derivatives with respect to coordinate x. These arc spin fluxes 
of mobile electron and lattice spins. After integration in @, they do not de- 
pend on Xf and remain finite in the limit — > 0. Finally, the third summand 
appears, which is proportional to a singular ^-function inside the SB layer. This 
summand arises due to nonequilibrium longitudinal spin injection by current. 
After integration in 0, it does not depend on Af also. 

Going the way indicated, we obtain the following vector boundary condition: 

J{X F + e) - J{-e) + J m {\f + e) - J M (-e) + (G sd ) =0. (7) 

The mobile electron spin flux J in the condition J2J originates from the 
equation @. The lattice magnetization flux 3m follows from the intralattice 
exchange effective field of Eq. Q . The following expression being valid: 

J m = 4M,S (8) 
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and a = yAM is a lattice magnetization diffusion constant. 

To calculate the last term (G s d) in Q, we take the sum equation for M + m 
and get 



(G sd ) = j([M, H sd ]) + ya([m, M]> = 7 a( M(.t 



(SAm/SM) dx' 



(9) 



where formulae and Q are used under the calculations. The condition J7J 
remains valid at x = L, if we replace: Xp + e — ► L + e, — e — > L — e and put 
(G sd ) -> 0. 



5 Nonequilibrium mobile electron spin flux and 
magnetization 

To make the system of equations and boundary condition be well 

defined, we should express all the quantities involved via the vectors M and m. 
As it may be seen from (J7J, we should consider the spin flux J outside the SB 
layer only. In the region x > Xp electrons occupy subbands having spins parallel 
to M (|) and antiparallcl to M (J,). Therefore, the following representations 
are valid: m = /is(n-f — rij)M = mM and J = (/is/e)(j| — j'jJM, where fig is 
Bohr magneton, n\ ^ and are partial electron densities and current densities 
in the spin subbands, respectively. The total electron density n = rif + and 
j — jf + do not depend on x and t because of local neutrality conditions 
in metal and one-dimensional geometry of our model. We use the well-known 
"diffusion-drift" formula for j'f ^ to express the partial currents via j, n and m. 
The corresponding calculations are direct and were done completely in |S] . The 
final result for x > Xf may be written as follows 

where Q = (erf— may be understood as an equilibrium conductivity 

spin polarization parameter with a^i = Mf,i n T>l' an< ^ ^ ~ ( a ] ^l + cr l^T)/( (7 T + 
<Ti) is the effective spin diffusion constant, quantities /ij^ and -Df ,x being partial 
mobilities and diffusion constant, respectively. To obtain ()10(l . an additional 
assumption should be made [B], namely, j/jo <C 1, where jd = enl/r is a 
characteristic current density in the layer 2. With typical parameter values, 
n ~ 10 22 cm" 3 , i~3x 1CT 6 cm, r ~ 3 x lCT 13 s, we get j D ~ 1.6 x 10 10 
A/cm 2 . 

The lattice in the layer 1 is pinned. Therefore, the magnetization flux J(— e) 
may be written similar to (|l(Jf> but without the spatial derivative. We introduce 
Mi = Mi/Mi and have 

J(- £ ) = ^QijM 1 . (11) 

e 

The latter flux has longitudinal and transversal components, which are: 
J||(-e) = ^Q 1J (M 1 M( £ ))M(e) 
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and 



(12) 



Mb, 



J±(-e) = ^Q 1 j M(e), M lt M(e) 
e 



We should calculate now the noncquilibrium magnetization for x > Xf. In 
the region an effective frequency of the motion cu is determined by the precession 
in relatively small fields: u> ~ jH, jHd, jH a -C 7«M = uj s d- We assume the 
conditions wt < 1 and w S (jT ^> 1 are valid for typical values r~3x 1CT 13 
s, a ~ 2 x 10 4 , M ~ 10 3 G. It allows neglect the time derivative in JSJ and 
substitute the formula l|10|) for flux J. Then the equation reduces to 



dx 2 



Ar 



= 0, 



(13) 



where / = v Dt is the spin diffusion length. Analogous relations may be written, 
of course, for the nonmagnetic layer 3 with the following modifications: 777,3 — 
and Q3 = 0. 

The next step is to find the solution of (13) that satisfies the conditions: 1) 
J|| (— e) = J (Xf + e), 2) J (L — e) = J (L + e), 3) continuity of subband chem- 
ical potential difference at the interface point x = L (see ^1 El for details). 
This solution takes the form 



Am(x) 



J 



fjL B n 



jr> sinh A + v cosh A 
x [cosh(A - + ^sinh(A -£)]}> 



{Qcoshe+[Qi(M 1 M(e))-Q] 



(14) 



where A = L/l, £ = x/l and parameter v characterizes the influence of the layer 
3 (typically v ~ 1, see pHE] f° r details). We substitute (|T4l into expression 
(0 and use the variational derivative (5(M 1 M(e))/JM(x) = M^(i — e). Then 
we get 



(Gad) = la^BnQil — 
JD 



1 



sinh A + v cosh A 



M(e),M 1 



(15) 



because of (S(x — e)) = 1. 

Turning back to condition J7J) , we see all the terms are now explicitly depen- 
dent on the vectors M and m. The latter vectors satisfy the equation of motion 
Therefore, the way opens now for solving our problem directly. 



6 Initial steady state and fluctuation instability 

To illustrate the possibilities of solution, we assume simple situation where field 
H is applied along the positive direction of z axis, H a is parallel to this axis 
too and the angle x = 7r - I n the situation all the equations and boundary 
conditions are satisfied for the following steady state magnetization vector M: 
M x = My = 0, M z = 1, and Mi = -z. 



G 



We introduce the fluctuations AM by the equality: M = z + AM with 
MAM = and AM = |AM| < 1. We linearize the equations LLG and 
condition Q with respect to AM and lay = — 4nAIAM x x. Then we get 
the equations for fluctuations AM XtV ~ exp(iqx — iuot): 



d 2 AM x (O, — ikuj) A - iu) A - 

5 a ~ — -AM, AAf a = 0, 

aar a a 



(16) 



d 2 AM y (tl y - inuj) 



IjJ 



dx 2 a v a 



AM, H AM, = 0, 



where two characteristic frequencies appear: fi^ = j (H + H a +47rM), Q y = 
7 (iJ + H a ) and transversal part of the vector boundary condition {7J) for fluc- 
tuations at x — takes the form: 

?^h- = kAMy — pAM x , 

(17) 

9AM, 



= —kAM x - pAMy. 

Parametersfc = (^BjQi)/{e^M 2 A) andp = ol[Ib3Q\t\Y— ^/(sinh A+^cosh A)]/(eAM) 
are proportional to current density j and describe respectively the effects of 
spin-transfer torque and spin injection. 

Transversal terms of the boundary condition {7|) at x = L become: 

^ = 0, ^=0. (18) 
ox ox 

Solving the standard boundary problem 1161) <|18[) . we come to the charac- 
teristic wave number 

1 



q 2 = — [n x + n y - 2iKLu ± ^J(n x - n y ) 2 + 4w 2 ] (19) 

and dispersion relation 

qLt&nqL = —pL ± ikL (20) 
for determining the spin- wave spectrum and decrement co(q). 



7 Steady state arising due to instability devel- 
opment 

The dispersion relation Q2(Jf) is the same as in our previous works ^3 and 
therefore all the results about the instability thresholds and current dependent 
spectrum remain valid. But in our new approach, contrary to |13l I14j . we have 
much more simple equations of motion, which do not contain any singular terms. 
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It allows solve more complicated problems now. One of such problem is: what 
will be a steady state the instability development leads to? 

We confine ourselves to the case of large enough thickness L > I, \f~A cor- 
responding to the injection mechanism dominated. In such case, the only polar 
angle 0(x) exists that determines the stationary magnetization M orientation 
and the problem reduces to solving of nonlinear equation 



b ——, j — smftcosi 
dx z 



--|-sin0 = O, 

tla 



(21) 



where 5 = y/A/fl. Transversal vector boundary conditions in the polar 
coordinates are 



d9 

— — hp sin t 

dx 







x=0 



, dO 
and — 

dx 



0. 



x—L 



(22) 



Equation (|21|l describes a "classical" particle moving in a potential field. 
As it is well known, such type of equation may be integrated easily. At pS < 

1/2 

(1 + H/Ha) 1 , there is no instability and only one solution 6 — exists. Under 
the opposite condition, the instability occurs and a nonzero solution appears 
which describes the new steady state orientation angle 6(x) in the layer 2: 

0(x) = 2arccosJ ((hj/jo) t&nh(x/h5) + l)(l - h 2 ) 1 / 2 



(hj/jo + tanh(x/M)) 2 - h 2 ((hj/j Q ) t&nh(x/h5) + l) 



-1/2 



», (23) 



where h = (l + H/H a ) 



-1/2 



j a = eSH a /afi B TQi. 
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